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ABSTRACT 

We present a model for the evolution of the galaxy ultraviolet (UV) luminosity function (LF) across 
cosmic time where star formation is linked to the assembly of dark matter halos under the assumption 
of a mass dependent, but redshift independent, efficiency. We introduce a new self-consistent treatment 
of the halo star formation history, which allows us to make predictions at z > 10 (lookback time < 500 
Myr), when growth is rapid. With a calibration at a single redshift to set the stellar-to-halo mass 
ratio, and no further degrees of freedom, our model captures the evolution of the UV LF over all 
available observations (0<z;<10). The significant drop in luminosity density of currently detectable 
galaxies beyond z ~ 8 is explained by a shift of star formation toward less massive, fainter galaxies. 
Assuming that star formation proceeds down to atomic cooling halos, we derive a reionization optical 
depth T = 0.056ito oioj fully consistent with the latest Planck measurement, implying that the universe 
is fully reionized at z = T.Sd^p'gg. In addition, our model naturally produces smoothly rising star 
formation histories for galaxies with L < L* in agreement with observations and hydrodynamical 
simulations. Before the epoch of reionization at z > 10 we predict the LF to remain well-described 
by a Schechter function, but with an increasingly steep faint-end slope (a ~ —3.5 at z ~ 16). Finally, 
we construct forecasts for surveys with JWST and WFIRSTand predict that galaxies out to z ~ 14 
will be observed. Galaxies at z > 15 will likely be accessible to JWST and WFIRST only through 
the assistance of strong lensing magnification. 

Subject headings: cosmology: theory, galaxies: high-redshift, stars: formation 


1. INTRODUCTION 

The rest-frame ultraviolet (UV) galaxy luminosity 
function (LF) and its evolution with redshift are cru¬ 
cial tracers of galaxy properties over cosmic time. In 
particular, UV light can be used efficiently to measure 
the star formation rate (SFR), because photons at rest- 
frame wavelengths around 1500 A are primarily produced 
by young, massive, and short-lived stars. Current obser¬ 
vations characterize the UV LF over the large majority 
of the history of the universe, ranging from studies in 
the local universe from Galaxy Evolution Explorer data 
(e.g., Burgarella et al. 2006 to Hubble Space Telescope 
Wide Eield Camera 3 {HST /WFC3) observations which 
now reach redshift z ^ 8 —10, i.e. lookback times greater 
than 13 Gyr (e.g. Bouwens et al. 2015b). Transforma¬ 
tional results on the rest-frame UV light from the epoch 
of reionization, when the universe was less than 0.8 Gyr 
old, have been made possible by large and dedicated com¬ 
munity efforts that have identified a sample of more than 
1000 galaxy candidates at z > 7, spanning a large range 
in luminosities. This progress is thanks to a variety of 
surveys, including ultradeep observations in blank fields 
with the HUDF09 and HUDF12 campaigns (Bouwens 
et al. 2011; Dunlop et al. 2013; Illingworth et al. 2013), 
the use of cluster-scale lensing to probe intrinsically faint 
objects (e.g., Yue et al. 2014; Atek et al. 2015), the large 
area, panchromatic CANDELS survey (Grogin et al. 
2011; Koekemoer et al. 2011), and wide-held random¬ 
pointing surveys to identify more luminous but rarer ob- 
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jects at the bright end of the LF (e.g., Trenti et al. 2011; 
Bradley et al. 2012; Schmidt et al. 2014). 

The picture emerging from these observations is 
that the number density of galaxies decreases with 
increasing redshift, while the LF remains con¬ 
sistent with a Schechter (1976) form, 4)(L) = 

$*(L/L*)“ exp {—L/L*)/L* out to z ~ 8, albeit the un¬ 
certainty and co-variance in the best htting parameters 
grows at high redshift (Schmidt et al. 2014; Bouwens 
et al. 2015b). However, other studies hnd either an 
equally good ht with a double-power law at z ~ 8 
(Finkelstein et al. 2015b), or a preference for the dou¬ 
ble power law at z ~ 7 (Bowler et al. 2014; Bowler et al. 
2015; but note that the z ^ 7 investigations are based 
on ground, rather than space observations). Irrespective 
of the form of the LF at the bright-end, whose evolution 
might be linked to changing astrophysical conditions of 
high-redshift galaxies, such as reduced feedback at early 
times (Somerville et al. 2008), it is clear that overall, 
the observed population of galaxies {M^b ^ — 17) at 
z > 7 does not produce sufficient photons to ionize the 
universe. This is consistent both with theoretical mod¬ 
eling (Trenti et al. 2010; Raicevic et al. 2011; Robertson 
et al. 2015) and indirect probes, such as the lack of detec¬ 
tions of gamma ray burst host galaxies (Trenti et al. 2012; 
Tanvir et al. 2012), which suggest that even the deepest 
HST observations are seeing only the tip of the iceberg of 
the population of star-forming galaxies. Therefore, it is 
likely that a very faint population of unseen dwarf galax¬ 
ies at z > 6 is the main contributor to the UV luminosity 
density and in turn to the ionizing photon budget (e.g.. 
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see Alvarez et al. 2012). Indirect observational support 
for extrapolation of the UV LF is provided by the de¬ 
tection of ultrafaint galaxies behind gravitational lenses 
at z ~ 2 (Alavi et al. 2014). However, see Giallongo 
et al. (2015) and Madau & Haardt (2015) for a discus¬ 
sion of the potential contribution of active galactic nuclei 
to reionization. 

HST is both wavelength and aperture limited to obser¬ 
vations at z < 11, but future progress to characterize the 
properties of the UV LF at high redshift will be boosted 
by the upcoming James Webb Space Telescope {JWST; 
Gardner et al. 2006) and Wide-Field Infrared Survey 
Telescope ( WFIRST; Spergel et al. 2015). These ob¬ 
servatories are expected to extend the frontier of galaxy 
detection to before the epoch of reionization at z > 10, 
when the first generation of galaxies were being assem¬ 
bled and hydrogen in the universe was predominantly 
neutral. 

Motivated by this upcoming improvement in discov¬ 
ery capabilities of high redshift galaxies, we aim here 
at predicting the UV LFs at z > 10. For this we use 
a simple, yet successful, semi-analytic framework intro¬ 
duced previously by Trenti et al. (2010, 2015); Tacchella 
et al. (2013) which assumes that the main driver of the 
evolution of the galaxy UV LF is the growth and hier¬ 
archical assembly of dark matter (DM) halos. Whilst 
this simple approach does not have the power of numer¬ 
ical hydrodyamic simulations (e.g., Genel et al. 2014; 
Furlong et al. 2015) to make predictions about galaxy 
properties such as morphology, and the role played by 
satellite galaxies, it provides a simple, robust and empir¬ 
ically calibrated method to make predictions for the evo¬ 
lution of global galaxy properties. In particular, the sim¬ 
plicity of our approach avoids the degeneracies of large 
multi-parameter numerical simulations, and allows us to 
calibrate at one reference redshift the complex physics 
that regulates the conversion of baryons into stars, and 
then to focus on investigating how DM halo assembly 
drives redshift evolution under the assumption that star 
formation efficiency is redshift independent at fixed halo 
mass. 

The key assumptions of the framework are that, for 
star forming galaxies, the stellar mass content of a DM 
halo depends on the halo mass but not on redshift, 
and that the stellar mass built up has a characteristic 
timescale given by the halo assembly time. The mass 
dependent efficiency of converting halo mass into stellar 
mass is calibrated empirically at one redshift, where good 
observational constraints are available (e.g. z = 4 — 5), 
and then applied to predict the redshift evolution of the 
UV LF by combining evolution in the DM halo mass 
function (HMF) and halo assembly time with stellar pop¬ 
ulation synthesis modeling. Since halos assembled more 
rapidly in the past, the SFR was higher at high red¬ 
shift, explaining the general brightening of the luminosity 
versus halo mass relation found by abundance matching 
studies (e.g. Gooray & Milosavljevic 2005). 

This simple strategy has been remarkably successful 
in describing the UV LF evolution with redshift (Tac¬ 
chella et al. 2013), with results similar to those obtained 
by other studies based on matching DM halos to galaxy 
luminosity (e.g. Lacey et al. 2011; Munoz 2012; Jaacks 
et al. 2013; Behroozi & Silk 2014; Dayal et al. 2014; 
Mashian et al. 2015). However, the earlier implemen¬ 


tations had two limitations. First, the model took into 
account only the star formation happening in the second 
half of the halo assembly history: the time the halo took 
to grow from Mh/2 to Mh- In addition, the calibration of 
the star formation efficiency was not guaranteed to be in¬ 
ternally self-consistent with the halo assembly history of 
a DM halo over redshift. Given that the UV luminosity 
of a galaxy is only weakly sensitive to the star forma¬ 
tion history of stellar populations with ages greater than 
a few hundred Myr (Madau & Dickinson 2014), these 
limitations had very little impact at z < 8. However the 
situation is potentially different at early times, when the 
assembly time goes below 100 Myr. Thus, the previous 
version of the model could not be trusted to formulate 
predictions at z > 10. 

In this paper, we develop and present an improved and 
self-consistent model that describes the full star forma¬ 
tion history of a DM halo, and we verify that (1) it con¬ 
tinues to describe well the UV LF at z < 8; (2) it success¬ 
fully reproduces the latest observations at z > 8, provid¬ 
ing an important validation of the approach introduced 
before such observations were available. In addition, we 
make detailed predictions for future surveys at z > 10 
with the James Webb Space Telescope (JWST) and the 
Wide-Field Infrared Survey Telescope (WFIRST), and 
we discuss the implications of our results in the context 
of the ionizing photon budget. 

This paper is organized as follows: in Section 2 we 
introduce our new model and its calibration; Section 3 
describes our model results, and our predictions for fu¬ 
ture surveys; and we summarize and conclude in Sec¬ 
tion 4. All magnitudes are AB magnitudes and we use 
the Planck Gollaboration et al. (2015) cosmology, with 
Dm = 0.315, Da = 0.685, Db = 0.0490, h = 0.6731, 
as = 0.829 and n* = 0.9655. 

2. MODEL DESCRIPTION 

Our model considers the growth of DM halos to be 
the most important driving force in the growth of galax¬ 
ies, and aims at predicting the evolution of the UV lu¬ 
minosity function with a minimal number of assump¬ 
tions. We thus assume that the SFR is proportional 
to: (1) the halo mass, through a mass-dependent but 
redshift-independent efficiency £(M/j), which is the ra¬ 
tio of the stellar mass formed during the halo assembly 
time ta{Mfi,z) to the final halo mass; and (2) to the in¬ 
verse of the halo assembly time. The halo assembly time, 
ta(Mh, z), of a halo of mass Mh observed at redshift z is 
the lookback time at which the progenitor halo had mass 
Mh/2 (Lacey & Gole 1993), and decreases at higher red¬ 
shift. This implies that DM halos at fixed mass host star 
forming galaxies with stellar mass which is independent 
of redshift, but with stellar populations that are younger 
at higher redshift. 

Since the HMFs and assembly times are fully defined 
by the cosmological model parameters, these assump¬ 
tions allow us to calibrate our model at one redshift 
to derive e{Mh), and then construct predictions for the 
galaxy luminosity function at all other redshifts from the 
DM HMF and the halo assembly time. 

2.1. Star formation prescription 

The UV luminosity of a galaxy most strongly depends 
on its youngest stars, while stellar populations older than 
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a few hundred Myr contribute little. However, at 2 : <8, 
the halo assembly time is less than 100 Myr, so to predict 
the UV luminosity at high redshift it is necessary to con¬ 
sider periods of star formation before the halo assembly 
time. 

To include multiple epochs of star formation, as a halo 
grows in DM mass, we define the star formation history 
of a halo as a linear combination of constant bursts nor¬ 
malized by the length of each burst. Thus we define the 
SFR for a halo with mass Mh in each epoch between 
times ti and as: 


SFR{ti, U+i,Mh) = MhX 


sjMtm 
2*(t*+i - U) 


( 1 ) 


where we define to as the lookback time for a halo ob¬ 
served at redshift 2obs and ti>o = ta{Mh/2'‘~^, Zi-i), 
where to {za) is the halo assembly time (redshift). We 
similarly define zq = Zohs and 2i>o = Zi_i). 

We calculate the halo assembly time as defined by Lacey 
& Cole (1993) in the extended Press-Schechter formal¬ 
ism (Bond et al. 1991) using an ellipsoidal collapse 
model (Sheth et al. 2001; Giocoli et al. 2007). We use the 
median of the probability distribution of assembly times 
for each halo. While this assumption does not take into 
account variations in the luminosity of individual galax¬ 
ies, there is a minimal effect on the global LF from ne¬ 
glecting scatter in halo assembly times, as demonstrated 
by Tacchella et al. (2013). 

We define the redshift-independent efficiency of star 
formation, e{Mh), as the ratio of the stellar mass formed 
during the halo assembly time to the final halo mass. 
Thus, to make predictions, we only need to calibrate 
e[Mh) at one redshift (see Section 2.3 and Figure 1) and 
can use the derived e{Mh) for all further predictions. 

Figure 2 shows the star formation history of halos of 
fixed final mass Mh = 10^^ Mq observed at 20 = 2, 6 
and 10, calculated using the SFR in Equation (1). As 
e{Mh) is redshift independent, these halos will also have 
identical stellar masses at their observed redshifts. The 
SFR shown in Figure 2 increases in each epoch as the halo 
grows from Mh/M to Mh, because e{Mh) decreases with 
decreasing halo mass (for Mh< 10^^ -M©) more rapidly 
than the shortening of the halo assembly times as the 
lookback time grows. This behavior of our model is fully 
consistent with strong evidence of rising star formation 
histories with redshift from both numerical simulations 
and observations (Finlator et al. 2011; Papovich et al. 
2011; Jaacks et al. 2012; Behroozi et al. 2013b; Lee et al. 
2014). Thus the greatest contribution to the stellar mass 
is during the halo assembly time as the halo grows from 
Mh/2 to Mh- This figure also illustrates how the short 
halo assembly times at high redshift require a consider¬ 
ably higher SFR to form the same final stellar mass. 

We include the contribution from star formation in suc¬ 
cessively smaller halo progenitors by summing the terms 
from Equation (1). The sum is truncated when the pro¬ 
genitor halo mass is below the cooling threshold, i.e. 
e{Mh < 10®Mq) = 0. 

Thus we can derive the stellar mass as: 


M^{Mh) =MhX 


g(Mjy) 

2* 

i=0 


( 2 ) 



Fig. 1.— The efficiency of star formation, the ratio of stellar 
mass formed during the halo assembly time to halo mass, see Equa¬ 
tion (2), derived at the calibration redshift z ~ 5, as described in 
Section 2.3. The shaded region shows Icr confidence range. 


which is redshift independent. 

To compute the UV luminosity of a halo we popu¬ 
late every halo with a galaxy with a stellar population 
based on the simple stellar population (SSP) models of 
Bruzual & Chariot (2003). We assume a Salpeter ini¬ 
tial mass function (IMF) between O.IM© and lOOM©, 
as low mass stars do not contribute much to UV lumi¬ 
nosity, and constant stellar metallicity Z — 0.02Z©. We 
neglect redshift evolution in metallicity as the UV lu¬ 
minosity does not depend strongly on metallicity under 
the assumption that current and future HST/JWST sur¬ 
veys detect galaxies living in relatively massive DM halos 
{Mh > 10® M©) where multiple generations of early star 
formation began enriching the gas at 2 > 20 —40 (Bromm 
et al. 2009; Trent! & Stiavelli 2009; Smith et al. 2015), so 
that by 2 < 16 the metallicity has risen to Z>0.01 Zq. 
We define £bc( 0 as the luminosity at 1500 A of an SSP 
of mass IM© and age t. The total UV luminosity of a 
galaxy observed at redshift z is obtained by integrating 
over the SFR (Equation (1)) and SSP luminosity in each 
epoch of star formation: 

^ ^ pti+1 

L{Mh,z)=Y. SFR{U,U+i,Mh)e^c{t)dt (3) 

z=0 

Where U are the halo assembly times for the half-mass 
progenitors defined above. 

2.2. Dust extinction 

The observed UV luminosity is significantly attenua¬ 
tion by dust extinction, particularly at 2 < 4. Thus, we 
include dust extinction in our model, following closely 
the procedure adopted in observations of Lyman-break 
galaxies. We assume a spectrum modeled as f\ ~ A^, 
and extinction Auv = 4.43 -I- 1.99/3 (Meurer et al. 1999). 
Following Trent! et al. (2015) and Tacchella et al. (2013) 
we model the observations of (3 by Bouwens et al. (2014) 
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Redshift, 2 : 



Fig. 2.— The star formation history, as described by Equation (1) 
in our model, of a halo of fixed mass 10^^ if observed at 

20 = 2, 20 = 6 or 20 = 10. We label the halo assembly time for 
the 20 = 2 halo and the mass of the halo progenitor at the start of 
each constant star formation epoch. 


as: 


(/?(Z, Myv)) — 

i/^Moiz) - c)exp 


g^(z)[MoV-Mo] 

Pmo{z)-c 


— Mq] + Pmo (z) 


+ C 


(4) 

Muv > Mq 
Muv < Mq 


where c = —2.33, Mq = —19.5, and the values of Pmq 
and d/3/dMQ are taken from Table 3 in Bonwens et al. 
(2014) and linearly extrapolated to lower and higher red- 
shifts. Using this linear plus exponential model for jd we 
fit the Bouwens et al. (2014) observations well - which 
show evidence for a curved relation between fd and Muv 
The exponential fit at faint magnitudes avoids unphysical 
negative dust corrections. 

We assume a Gaussian distribution for jd at each Muv 
value (with dispersion ap = 0.34), giving the average 
extinction (z4uv) = 4.43-1-0.79ln(10)(T^-|-1.99(/l). We use 
this average extinction to calculate observed UV fluxes 
in our model. 


2.3. Calibration 

We calibrate the model by finding e{Mh) such that 
L{Mfi,Zc) = Zc), where Zc is the calibration 

redshift. To find Zc) we derive an empirical 

relation between observed luminosity and halo mass by 
performing abundance matching (Mo et al. 1999) be¬ 
tween the HMF, n{Mh,z), and observed LF, $(L, z), 
assuming every halo hosts one galaxy: 

pOO pOO 

/ n{Ml,z)dMl= ^L',z)dL' (5) 

where we use the Sheth & Tormen (1999) HMF and cal¬ 
ibrate at z ~ 5 using the Bouwens et al. (2015b) LF. 
We calibrate over the halo mass range 10^ — lO^^M©, 


and therefore extrapolate the LF over Muv = —7.5 to 
Muv = —24.5 to perform the abundance matching. 

We can then calculate the UV luminosity as a function 
of halo mass at any redshift using Equation (3). The 
largest contribution to the uncertainty in the predicted 
luminosity is from the uncertainty in the calibration LF; 
we plot 1(7 confidence bounds from this uncertainty in 
the plots which follow. 

3. RESULTS 

3.1. Stellar masses and ages 

Figure 1 shows the efficiency of star formation as a 
function of halo mass, e{Mh), derived at the calibration 
redshift z ^ 5. Halos with mass 10^^ — 10^^ Mq 

have the highest star formation efficiencies, which is con¬ 
sistent with the consensus picture that low-mass halos 
have lower efficiencies because of supernova feedback, 
while star formation in higher mass halos is affected by 
strong negative AGN feedback. Our derived e{Mh) is 
consistent, to first approximation, with the results ob¬ 
tained by Behroozi et al. (2013a) who calculated the ra¬ 
tio of stellar to halo mass at different redshifts (z < 4) 
by applying abundance matching techniques between the 
stellar mass function and HMF. We note that e{Mh) is 
redshift-independent in our model, but that the evolution 
of the HMF and halo assembly times allows to us predict 
the UV luminosity using Equation (3), via our rising star 
formation history (see Figure 2). Thus we find that with 
our simple redshift-independent e{Mh) the star forma¬ 
tion at high redshift proceeds much more rapidly than 
at low redshift. 

With Equation 1, we can derive the model predictions 
for average galaxy ages as a function of halo mass and 
redshift, which are shown in Figure 3. Galaxy ages are 
challenging to constrain observationally, due to the de¬ 
generacy between dust extinction, age and metallicity in 
spectral energy distribution (SED) fitting, which are es¬ 
pecially severe at z>4 where HST covers only the rest- 
frame UV wavelengths. Our modeling results are con¬ 
sistent with the picture emerging from multiple recent 
studies that combine HST and Spitzer/IRAG data and 
find that the majority of z > 4 galaxies to have old stel¬ 
lar populations (> 100 Myr) and relatively low specific 
SFR (~ 100 Mq yr“^) (Oesch et al. 2013a; Straatman 
et al. 2014; see also Gonzalez et al. 2011 for earlier stud¬ 
ies reaching the same conclusion but without account¬ 
ing for nebular emission lines). However, other inde¬ 
pendent studies based mostly on overlapping datasets 
reached the different conclusion that z>4 galaxies have 
a high chance of being young systems (ages < 50 Myr) 
with high specific SFRs (de Barros et al. 2014; Finkel- 
stein et al. 2015a), suggesting that improved observa¬ 
tional constraints are needed to evaluate the hdelity of 
stellar ages predicted by our model. 

The contributions to the total UV luminosity at fixed 
halo mass from two epochs of star formation as a func¬ 
tion of observed redshift are shown in Figure 4. We find 
that the earlier epoch of star formation, as the halo grows 
from Mh/4 to Mh/2 adds a negligible contribution to the 
total UV luminosity, suggesting that recent star forma¬ 
tion is the most important contribution to the UV lumi¬ 
nosity even at z > 10 when the assembly time is short. 
This is consistent with recent clustering studies (Barone- 
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Fig. 3. — Average stellar population age in our model and halo 
assembly time (Lacey & Cole 1993) as a function of halo mass and 
observed redshift. 



Fig. 4.— Contribution to predicted UV luminosity from the 0th 
and 1st order terms in Equation (3) as a function of observed 
redshift for fixed halo mass, = 10^^Mq. The contribution 
from star formation during the halo assembly time (solid) domi¬ 
nates, with the contribution from the earlier star formation epoch 
(dashed) increasing with redshift. The SFR is constant in both 
epochs. 

Nugent et al. 2014) which found a high duty cycle for 
galaxies at z > 6.5, possibly due to their bright, young 
stellar populations, and with the prediction of smoothly 
rising SFRs from numerical simulations (Finlator et al. 
2011) and observations (Papovich et al. 2011; Behroozi 
et al. 2013b; Lee et al. 2014). Thus, we only consider the 
contribution from the first two terms in Equation (3) in 
predicting the UV luminosity. 

The stellar mass density as a function of redshift, ob¬ 
tained by integrating our model stellar mass functions 
(using Equation (2) to derive stellar masses) to a stellar 
mass limit of M* > 10 ®Mq is shown in Figure 5. We find 
good agreement with observations (Perez-Gonzalez et al. 
2008; Stark et al. 2013; Tomczak et al. 2014; Grazian 
et al. 2015; Song et al. 2015) taking into account the 
scatter in observations at z < 2, and the uncertain but 
potentially significant contribution of massive quiescent 
galaxies at low redshift. 


Fig. 5.— Stellar mass density (p*) obtained by integrating the 
stellar mass function derived from our model using Equation (2) 
to a stellar mass limit of M* > 10® Mq. We plot observations 
from Perez-Gonzalez et al. (2008); Stark et al. (2013); Tomczak 
et al. (2014); Grazian et al. (2015); Song et al. (2015). Shaded 
regions show the la confidence range. 

3.2. Luminosity functions and SFR density 

The predicted observed L{Mh,z) for a range of red- 
shifts is plotted in Figure 6. Our model naturally pro¬ 
vides redshift evolution of L{Mh,z) through the evolu¬ 
tion of the HMF and halo assembly times. Dust extinc¬ 
tion (Section 2.2) significantly affects the high mass end 
of the relation at low redshift. At z > 10 we see the 
high mass end does not evolve much with redshift, moti¬ 
vating the model of Mashian et al. (2015) which uses an 
empirical redshift-independent L{Mh). However, there is 
significant evolution at lower mass, which comprises the 
greatest contribution to the photon budget available to 
reionization because of the steep faint-end slope. 

The model UV LFs at z < 7 and z > 7 are shown 
in Figures 7 and 8 respectively. The model is remark¬ 
ably consistent with the observed data (Arnouts et al. 
2005; Oesch et al. 2010; Oesch et al. 2013b, 2014; Alavi 
et al. 2014; Bowler et al. 2015; Finkelstein et al. 2015b; 
Bouwens et al. 2015b,a), which is expected due to the 
success of our previous implementation of this class of 
models (Tacchella et al. 2013; Trent! et al. 2010, 2015). 
We find the model marginally overpredicts the bright end 
of the LF at z ^ 0.3 and z ^ 2 by ^ 0.1 dex due to 
the difficulty in modeling dust extinction at these red- 
shifts, but the observations are still within 2cr of our 
model. In particular, our model predicts a steepening 
of the faint-end slope at higher redshifts, consistent with 
the observed trend (Bouwens et al. 2015b). 

At z > 7, the model describes the most recent ob¬ 
served data (Oesch et al. 2013b, 2014; Finkelstein et al. 
2015b; Bouwens et al. 2015b,a) well, validating our sim¬ 
ple approach. Our model predicts the trend of steepening 
faint-end slope to continue at z > 10, and number den¬ 
sities to drop rapidly. With JWST capabilities, except 
in an extremely wide-held survey, the UV LF at z > 10 
will be observed as a steep power-law function. This is 
agreement with the semi-analytic results of Behroozi & 
Silk (2014) who hnd a steepening power law slope at faint 
magnitudes and a signihcant drop in number densities. 

For comparison, we also plot in Figure 9 the model LFs 
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TABLE 1 

Best-fit Schechter parameters for model LFs 


Redshift 

OL 

M* 

log(4>*[mag ^Mpc ^]) 

2 ~ 0 

-1.68 ±0.09 

-19.9 ±0.1 

^•^'+0.08 

2 ~ 2 

-1.46 ±0.09 

-20.3 ±0.1 

9 

^•^^+0.09 

2 ~ 4 

-1.64± 0.11 

-21.2 ±0.2 

— 2 

2 ~ 5 

-1.75 ±0.13 

-21.2 ±0.2 

o -| Q — 0.15 

— 24 

2-^6 

-1.83 ±0.15 

-20.9 ±0.2 

— S 

'^•-^^+0.25 

2 ~ 7 

-1.95 ±0.17 

-21.0 ±0.2 


2 ~ 8 

-2.10± 0.20 

-21.3 ±0.4 

—4 

72 

2 ~ 9 

-2.26 ±0.22 

-21.2 ±0.4 

—4 e:n—0-29 
^•^^+1.36 

2 ~ 10 

-2.47 ± 0.26 

-21.1 ± 0.5 

-5.12 ±0.34 

z ~ 12 

-2.74 ± 0.30 

-21.0 ±0.5 

-5.94 ± 0.38 

2 ~ 14 

-3.11 ± 0.38 

-20.9 ±0.5 

-7.05 ±0.45 

2 ~ 16 

-3.51 ± 0.46 

-20.7 ± 0.6 

-8.25 ±0.51 


Note. - Fit performed between Mab = —17.5 and Mab = —22.5 


obtained by using the Finkelstein et al. (2015b) z ~ 5 LF 
for the calibration. The Finkelstein et al. (2015b) z ~ 5 
LF has a lower value of M* than that of Bouwens et al. 
(2015b), so it is not surprising that this calibration pro¬ 
duces slightly lower number densities compared to our 
reference using the Bouwens et al. (2015b) LF. How¬ 
ever, Figure 9 clearly shows that the overall evolution 
trends are unchanged and that the two calibrations pro¬ 
duce LFs which are consistent within one standard de¬ 
viation both between themselves and the observed data 
over 0 < z < 10. 

The best-fit Schechter (1976) function parameters for 
our LFs are shown in Table 1. The best-fit parameters 
are in good agreement with observations (Schmidt et al. 
2014; Bowler et al. 2015; Oesch et al. 2010; Finkelstein 
et al. 2015b; Bouwens et al. 2015b) given the large degen¬ 
eracies in Schechter function parameters. Encouragingly, 
we find the evolution of the derived Schechter parame¬ 
ters is in excellent agreement with the observed evolu¬ 
tion (Bouwens et al. 2015b; Bowler et al. 2015): we find 
dajdz ^ —0.1, dM*/dz ~ 0.1, and d\og{^*)/dz ^ —0.3 
between z ~ 4 and z ~ 8. We find the evolution of a 
and between z ~ 8 and z ^ 16 to be more dramatic: 
da/dz ~ —0.2, dM*/dz ~ 0.1, and d\og{^*)/dz ~ —0.5, 
consistent with the rapid evolution of ^ 10^° Mq halos 
in the DM HMF at these redshifts. 

Figure 10 shows the luminosity density and cosmic 
SFR density as a function of redshift. We calculate the 
luminosity density by integrating our model LFs down 
to a magnitude limit. We choose two hducial limits of 
-Miim = —17 (just fainter than current observational lim¬ 
its) and Miim = —12 (the theoretical mass limit for ha¬ 
los to cool). We calculate the SFR density, p* using 
the empirical relation from Madau et al. (1998) where 
SFR[MQ/yv\ = 8.0 x lO^^Ljergs/s/Hz] at 1500 A. We 
plot the densities and observations from Bouwens et al. 
(2015b) both with and without dust correction. At z < 8 
the observations are consistent with both magnitude lim¬ 
its, however the observations at z ^ 10 suggests a sig¬ 
nificant steepening of the relation at high redshift, as 
do results from numerical simulations (Genel et al. 2014, 
though with large uncertainty), which is consistent with 
our model with M\i^ = —17. The sample at z ^ 10 is 
limited however; more observational data at z > 8 are 
needed to confirm this result. 



Fig. 6.— The relationship between observed galaxy luminosity 
and halo mass as a function of redshift, plotted at 2 = 

2, 6, 8, and 16. Shaded regions show the Icr confidence range. 

3 . 3 . Forecasts for JWST and WFIRST 

We use our model to make forecasts for a representative 
set of JWST NIRCAM high-redshift dropout surveys us¬ 
ing the 5 near-IR filters. The surveys (properties sum¬ 
marized in Table 2) include an ultra-deep (UD) survey 
of 4 pointings (^ 40 arcmin^) exposed in 200 hours per 
pointing; a medium-deep (MD) survey of 40 pointings 
exposed in 20 hours per pointing; and a wide-held (WF) 
survey of 400 pointings exposed in 2 hours per pointing. 
We assume that the surveys will split the observing time 
so as to reach equal depth in all hve hlters, and estimate 
the limiting magnitude for an 8cr detection (in a single 
hlter) using the JWST Exposure Time Calculator. We 
also include the effects of gravitational leasing magni- 
hcation bias from strong leasing in blank helds, which 
is expected to distort the brightest end of high-redshift 
LFs (Mason et al. 2015; Wyithe et al. 2011). 

In Figure 11 we plot the predicted cumulative number 
counts for redshifts 8 < z < 16 and the regions acces¬ 
sible to these mock JWST surveys, as well as the re¬ 
gion accessible to WFIRST High-Latitude Survey (HLS, 
Spergel et al. 2015). The estimated number of dropouts 
are given in Table 2. 

Our model predicts a signihcant drop in number den¬ 
sity from z ~ 8toz ~ 10 compared to lower red¬ 
shifts (which is also seen in the observations, Bouwens 
et al. 2015b,a). The drop continues to high redshift, thus 
we find that no z ~ 16 galaxies would be detected in our 
mock JWST surveys. To detect 1 galaxy at z ~ 16 
in our UD survey would require ^ 40 pointings (^ 400 
arcmin^). We find that magnification bias in blank fields 
does not significantly affect our model even at the bright¬ 
est observable magnitudes at z > 10. The magnihcation 
bias effect is only noticeable in the exponential part of 
the LF, which is within reach only at z 8, but too weak 
otherwise for power laws with slope in the range -2 to -3.5 
(it is exactly neutral for faint end slope a = —2). Mason 
et al. (2015) showed the lensing effect was most signif¬ 
icant for a Schechter function LF at high redshift (see 
also Barone-Nugent et al. 2015). Thus we expect that 
without significant strong lensing, i.e. using galaxy clus- 
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Muv 

Fig. 7.— Predicted UV LFs at low (upper) and intermediate 
(lower) redshift. We show the LFs using the calibration (see Sec¬ 
tion 2.3) at 2 ~ 5 from Bouwens et al. (2015b), with Planck 2015 
cosmology (Planck Collaboration et al. 2015), with Planck 2015 
cosmology (Planck Collaboration et al. 2015). Points show the 
binned UV LFs and upper limits from Arnouts et al. (2005); Alavi 
et al. (2014); Oesch et al. (2010); Bouwens et al. (2015b); Finkel- 
stein et al. (2015b); Bowler et al. (2015); Atek et al. (2015). We 
note that the data from Atek et al. (2015) was made public af¬ 
ter our model was submitted and illustrates the consistency of our 
model with observations even at very low luminosity. Shaded re¬ 
gions show the 1(7 confidence range. 


ters as cosmic telescopes (e.g. the Hubble Frontier Fields 
Yue et al. 2014; Ishigaki et al. 2015; Coe et al. 2015; Atek 
et al. 2015), z > 15 is beyond the reach of JWST. 



Fig. 8.— Predicted UV LFs at high redshift. We show the LFs 
using the calibration (see Section 2.3) at z ~ 5 from Bouwens 
et al. (2015b), with Planck 2015 cosmology (Planck Collaboration 
et al. 2015). Points show the binned UV and upper limits LFs 
from Oesch et al. (2013b, 2014); Finkelstein et al. (2015b); Bouwens 
et al. (2015b,a). Shaded regions show the Icr confidence range. 



Fig. 9. — Predicted LFs at redshifts z ~ 2, 5, 10, 16 obtained by 
calibrating (see Section 2.3) our model with the Finkelstein et al. 
(2015b) LF at z ~ 5 (F15, dashed), compared to our reference 
calibration using the Bouwens et al. (2015b) LF at z ~ 5 (B15, 
solid). Shaded regions show the Icr confidence range, highlighting 
that within the uncertainty of the calibrations, the two approaches 
yield consistent results. 


3.4. Implications for reionization 

The timeline of cosmic reionization depends on the bal¬ 
ance between the recombination of free electrons with 
protons to form neutral hydrogen atoms, and the ion¬ 
ization of hydrogen atoms by Lyman continuum photons 
emitted by young stars. The UV luminosity density (and 
therefore, SFR density) at a given redshift allows us to 
calculate the number of photons available for reioniza¬ 
tion, and is most sensitive to the faint end of the LF. We 
can use this to infer the timeline of reionization by calcu¬ 
lating the ionized hydrogen fraction, Q{z)^ as a function 
of redshift given the following time-dependent differential 
equation: 


where hion is the comoving number density of ionizing 
photons, (iih) is the comoving number density of hydro¬ 
gen atoms, and the recombination time of the IGM (Sti- 
avelli et al. 2004; Robertson et al. 2015, and references 
therein) is: 

Fec{z) = [CaBiT)neil + zf] (7) 

where asiT) is the case B recombination (i.e. opaque 
IGM) coefficient for hydrogen, rig = (1 + Yp/AXp){nY{) 
is the comoving number density of electrons (assuming 
singly ionized He), Xp and Yp are the primordial hy¬ 
drogen and helium abundances respectively, and C = 
(n'jj) /is the “clumping factor” which accounts for 
inhomogeneity in the IGM. 

The production rate of ionizing photons can be related 
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TABLE 2 

Predicted Number Counts for Example JWST and WFIRST Surveys 


Redshift 

Dropout Filter 

UD (miim = 32.0) 

MD (miim = 30.6) 

WF (miin, = 29.3) 

WFIRST (miiAi = 26.5) 



~ 40 arcmin^ 

~ 400 arcmin^ 

~ 4000 arcmin^ 

2000 deg^ 

2 : ~ 8 

F115W 


548/11 

iOOO_5Q3 

61370/11111 

2 ~ 10 

F115W 

30/^7 

52/i 

102/51 

1026/505 

2 ~ 12 

F150W 

6/5 

10/5 

13/50 

47/15 

2 ~ 14 

F150W 

0.3/0/ 

0 

1 -h 

0 0 

0 

1 + 

0 0 

(M 

6 d 
+ 1 

d> 

2 ~ 16 

F200W 

0 

0 

0 

0 


Note. - Limiting magnitudes for a 8(j detection estimated with the JWST Exposure Time Calculator and WFIRSTHLS. 
The mock surveys are described in Section 3.3. These estimates include the boost from gravitational lensing 
magnification bias in blank fields (Mason et al. 2015; Wyithe et al. 2011). 




Fig. 10.— Luminosity density {pl) cosmic SFR density (p*) 
as functions of redshift, derived by integrating the model UV LFs to 
magnitude limits of Mab = —17 (green lines) and Mab = —12 (pur¬ 
ple lines). The dust corrected SFR densities for the two magnitude 
limits are shown as solid lines, dust uncorrected SFR densities are 
shown as dashed lines. The observed SFR densities from Bouwens 
et al. (2015b) are shown in black (dust corrected) and grey (dust 
uncorrected). Shaded regions show the Icr confidence range. 

to the total UV luminosity density, pl as 

hion — fesc^ionPL (S) 

where /esc is the average fraction of photons which escape 
galaxies to affect the IGM, and ^ion is the rate of ioniz¬ 
ing photons per unit UV luminosity, with units Hz/ergs, 
which depends on the initial mass function, metallicity, 
age and dust content of the stellar populations. There 
is an equivalent relation between hio„ and SFR den¬ 
sity (Madau et al. 1999; Shull et al. 2012), which requires 
the same stellar population modeling. 

All of the parameters involved are difficult to estimate, 
and may evolve with redshift as reionization progresses 
and the IGM evolves (Fuiianetto & Oh 2005; Shull et al. 
2012). In this work, we follow Schmidt et al. (2014) and 
use a distribution of parameters. For two limiting mag¬ 
nitudes (Mab = —17, corresponding to currently observ¬ 
able galaxies, and Mab = —12 , corresponding to atomic 
cooling halos) we assume the escape fraction is uniformly 
distributed between /esc = 0.1 — 0.3 (Ouchi et al. 2009), 
and we use a uniform distribution between (7=1 — 6 
for the clumping factor. Finally we model ^ion as a log¬ 
normal distribution with mean log^ion = 25.2 and stan¬ 
dard deviation 0.15 dex, which was the range obtained 
by Robertson et al. (2013) using the standard Bruzual 


Fig. 11. — Predicted number counts of galaxies brighter than 
apparent magnitude muv (rest-frame UV) per square degree for 
a range of redshifts based on our model LFs. We plot the cumu¬ 
lative number counts including the boost from gravitational lens¬ 
ing magnification bias (Mason et al. 2015; Wyithe et al. 2011) 
as solid lines, and without the magnification bias effect (dasheti 
lines). We plot the estimated coverage of future surveys as shaded 
regions: 3 mock JWST surveys detailed in Section 3.3 and the 
WFIRST High-Latitude Survey (Spergel et al. 2015). The calcu¬ 
lated number counts are given in Table 2. 

& Chariot (2003) models and measurements of the UV 
spectral slope by Dunlop et al. (2012). We assume an 
IGM temperature of 20, 000 K. 

Once the reionization history, Q{z)^ is known, an im¬ 
portant constraint is to compare the electron scattering 
optical depth with that inferred from CMB observations. 
The Planck Collaboration et al. (2015) reported a reion¬ 
ization value of T = 0.066±0.012, consistent with instan¬ 
taneous reionization at z = 8.8/]^'^. The optical depth 
as a function redshift is: 

T-iz) = aTne{l +z'fQ{z')j^^dz' (9) 

where c is the speed of light, ctt is the Thomson scatter¬ 
ing cross section and H{z) is the Hubble parameter. 

Figure 12 shows the reionization history: the ion¬ 
ized hydrogen fraction as a function of redshift, ob¬ 
tained by solving Equation (6) with our model luminos¬ 
ity density, sampling the distributions of input param¬ 
eters. Figure 13 shows the electron scattering optical 
depth as a function of redshift. For the LF magnitude 
limit Mab = —17, reionization is complete {Q = 1) by 
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Fig. 12.— The fraction of ionized hydrogen as a function of red- 
shift, obtained by solving Equation (6) with our model luminosity 
density. We plot our results from integrating the model UV LFs 
to two magnitude limits of Mae = —17 (green) and Mae = —12 
(purple), with Icr confidence regions as shaded regions. We also 
plot constraints derived from observations of: Lyo emission from 
galaxies (open circles, Ouchi et al. 2010; Pentericci et al. 2014; Tilvi 
et al. 2014; Faisst et al. 2014; Schenker et al. 2014); the Lya for¬ 
est (filled circles, Fan et al. 2006); the clustering of Lya emit¬ 
ting galaxies (square, Ouchi et al. 2010); GRB spectra damping 
wings (diamond, McQuinn et al. 2008); dark gaps in the Lya for¬ 
est (upper triangles, McGreer et al. 2015); quasar near zones (star, 
Venemans et al. 2015); and quasar spectra damping wings (lower 
triangle, Schroeder et al. 2013). We also plot the Planck Collabo¬ 
ration et al. (2015) redshift of instantaneous reionization. We note 
that the conversion from the Lya escape fraction to the global 
ionized hydrogen fraction is uncertain and relies on several model 
assumptions (Mesinger et al. 2015). 



Fig. 13. — The electron scattering optical depth, calculated using 
Equation (9) from our derived Q(z). We plot our results from 
integrating the model UV LFs to two magnitude limits of Mae = 
— 17 (green) and Mae = —12 (purple), with Icr confidence regions 
as shaded regions. We show the reionization optical depth value 
and its la confidence levels from Planck Collaboration et al. (2015) 
in grey. 

-^reion = 6.86lo|g, with T{zreion) = 0.0421°;°^. For the 
LF magnitude limit Mab = —12, reionization is complete 
by ^reion = 7-84to;98> with T(zreion) = 0.056to;o?o. 

The fainter magnitude limit, corresponding to atomic 
cooling halos of mass Mh ^ 10® Mq, is fully consis¬ 
tent with the Planck results, considering the uncertainty 
in the reionization model parameters. This calculation 


shows that ultrafaint galaxies can in principle provide 
enough photons to fully reionize the universe by z ~ 6 
to match observations of the Lya forest (Fan et al. 
2006). Both magnitude limits are broadly consistent 
with a range of constraints from observations, within 
the reionization model uncertainty: UV luminosity den¬ 
sities (Finkelstein et al. 2012) for observable galaxies; 
quasar near zones (Venemans et al. 2015); quasar spec¬ 
tra damping wings (Schroeder et al. 2013); GRB spec¬ 
tra damping wings (McQuinn et al. 2008); transmission 
(Fan et al. 2006) and dark gaps (McGreer et al. 2015) 
in the Lya forest; and the clustering of Lya emitting 
galaxies (Ouchi et al. 2010). 

Qualitatively, the non-negligible neutral fraction pre¬ 
dicted by our model at z>7 is consistent with the ob¬ 
served high optical depth of Lya (Ouchi et al. 2010; Treu 
et al. 2013; Pentericci et al. 2014; Schenker et al. 2014; 
Tilvi et al. 2014; Faisst et al. 2014, K. B. Schmidt et 
al. 2015, ApJ submitted), however the conversion from 
the Lya emission fraction to the volume filling factor 
of ionized hydrogen is difficult and requires several as¬ 
sumptions (Mesinger et al. 2015). In particular, to make 
constraints on reionization it is generally assumed that 
there are no changes in galaxy and the Lya emission line 
properties, which necessitates a rapid evolution of the 
global ionization fraction between z ~ 6 and z ~ 7. How¬ 
ever, recent studies have shown that the rapid decline in 
the Lya escape fraction at these redshifts cannot result 
only from the changing IGM attenuation (Mesinger et al. 
2015) but could also be explained by the co-evolution of 
the escape fraction of ionizing photons, /esc, (Dijkstra 
et al. 2014). Thus, the uncertainties in the ionization 
fraction from the Lya optical depth shown in our plot 
are likely underestimated, since they do not include these 
systematic effects. 

4. SUMMARY AND CONCLUSIONS 

We have presented a simple model for the evolution of 
the UV LF from 0 < z < 16, assuming that the average 
star formation history of galaxies is set by their halo mass 
and by the redshift (through the halo assembly time), so 
that halos of the same mass have the same stellar mass 
content independent of redshift. Our model builds upon 
previous similar implementations, but here we extended 
our framework to construct a self-consistent model which 
is capable of following the evolution of the star formation 
even when the halo assembly times become very short (at 

-^^lO)- 

Our key findings are as follow: 

1. Our model UV luminosity functions are very suc¬ 
cessful in matching observations at all redshifts 
where data are available (0<z<10). Overall, we 
find that the shape of the LF is well described by a 
Schechter function with faint-end slope increasing 
with redshift. This trend continues at higher red¬ 
shift, and we use the model to make predictions for 
LFs at z > 10, finding a faint-end slope a ~ —3.5 
at z = 16. 

2. Our model reproduces the observed cosmic SFR 
density well, indicating a sharp decline at z > 8 
with a magnitude limit of Mab = —17, consistent 
with observed data at z ~ 10. 
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3. Compared to previous more basic models (Tac- 
chella et al. 2013; Trenti et al. 2015), we find that 
the self-consistent inclusion of earlier periods of 
star formation does not significantly affect the total 
UV luminosity at a given halo mass and redshift, 
but it allows us to better reproduce the observed 
average stellar ages and stellar mass density of high 
redshift galaxies. 

4. Taking advantage of the ability of the model to 
make predictions at the earliest times, we inves¬ 
tigate the expected galaxy detections for future 
ultra-deep, medium-deep and wide-held surveys 
with JWST, and WFIRST. We predict that 2 ~ 14 
galaxies over a range of luminosities are in reach 
of these surveys. However, signihcant strong lens- 
ing magnihcation will be needed to push beyond 
z > 15. 

5. Finally, we investigate the implications of our 
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